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ABSTRACT 

We present numerical simulations aimed at exploring the effects of varying the sub-grid physics 
parameters on the evolution and the properties of the galaxy formed in a low-mass dark matter 
halo (~7x 10 10 fc. _1 M at redshift z = 0). The simulations are run within a cosmological setting 
with a nominal resolution of 218 pc comoving and are stopped at z = 0.43. For simulations that 
cannot resolve individual molecular clouds, we propose the criterion that the threshold density for 
star formation, usf, should be chosen such that the column density of the star- forming cells equals the 



threshold value for molecule formation, N 



M Q pc" 2 . In all of our simulations, 



an extended old/intermediate-age stellar halo and a more compact younger stellar disk are formed, 
and in most cases, the halo's specific angular momentum is slightly larger than that of the galaxy, and 
sensitive to the SF/feedback parameters. We found that a non negligible fraction of the halo stars 
are formed in situ in a spheroidal distribution. Changes in the sub-grid physics parameters affect 
significantly and in a complex way the evolution and properties of the galaxy: (i) Lower threshold 
densities tisf produce larger stellar effective radii R ei less peaked circular velocity curves V C (R), and 
greater amounts of low-density and hot gas in the disk mid-plane; (ii) When stellar feedback is modeled 
by temporarily switching off radiative cooling in the star forming regions, R e increases (by a factor 
of ~ 2 in our particular model), the circular velocity curve becomes flatter, and a complex multi 
phase gaseous disk structure develops; (iii) A more efficient local conversion of gas mass to stars, 
measured by a stellar particle mass distribution biased toward larger values, increases the strength 
of the feedback energy injection -driving outflows and inducing burstier SF histories; iv) If feedback 
is too strong, gas loss by galactic outflows — which are easier to produce in low-mass galaxies — 
interrupts SF, whose history becomes episodic; v) In all cases, the surface SF rate versus the gas 
surface density correlation is steeper than the Kennicutt law but in agreement with observations 
in low-surface brightness galaxies. The simulations exhibit two important shortcomings: the baryon 
fractions are higher, and the specific SF rates are much smaller, than observationally inferred values for 
rcdshifts ~ 0.4— 1. These shortcomings pose a major challenge to the SF/feedback physics commonly 
applied in the ACDM-based galaxy formation simulations. 

Subject headings: cosmology:dark matter galaxiesdormation galaxies:haloes 

methods:^— body simulations 
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1. INTRODUCTION 

The succesful development of cosmology in the last 
two decades, incarnated in the popular ACDM model, 
has provided a robust theoretical background for mod- 
eling galaxy f ormation and evo lution (for reviews see 
iBaugh l [20061: lAvila-Reesel [2007h . Due to the high 
non-linearity implied in the problem, cosmological N- 
body + hydrodynamical simulations offer the fairest 
way to attain suc h a modeling (for a recent review see 
iMaver et al.l l2008h . However, this method is hampered 
by the large -currently unreachable- dynamic range re- 
quired to model galaxy formation and evolution in the 
cosmological context (formally from the subpc scales 
of molecular cloud cores to Mpc scales, and from den- 



sities n 



10" 6 to > 10 2 



as well as by the 



complexity of the processes involved, mainly cooling 
and star formation (SF) and its feedback on the sur- 
rounding medium. Both of them occur at scales com- 
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monly well below the accesible resolution in simula- 
tions. Theoretical and observational results suggest that 
the star formation efficiency can be simulated with- 
out reaching the subparsec resolution if a simulation 
can resolve the mean density of giant molecu lar clouds 
(jGnedin et aLlEfflOl iKrumholz fc McKeell2005t) . unfortu- 
nately this seems not to be the case for SN and Stellar 
Wind feedback. Undoubtedly, modeling of SF and feed- 
back (subgrid physics) is of crucial importance for the 
simulations of galaxy formation and evolution. 

On the other hand, the SF rate (SFR) and its evolu- 
tion are among the most important observational prop- 
erties of galaxies, being basically responsible for coining 
the present-day galactic stellar populations. Star forma- 
tion is also crucial in terms of its feedback effect, which 
(i) regulates several thermo- and hydro-dynamical pro- 
cesses of the multi-phase ISM, including probably the SF 
itself; (ii) partially controls the galaxy assembly process 
through mass outflows and re-heating of the circumgalac- 
tic medium, and (iii) likely is also the driving agent of 
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the intergalactic chemical enrichment. 

In the context of the hierarchical ACDM scenario, the 
efficient radiative cooling of gas leads to the so-called 
"gas cooling catastrophe" : most of the available primor- 
dial gas is cooled, trapped, and transform ed into stars in- 
side t he smallest, early-collapsing halos (W hite fc Rees I 
119781 ). In order to avoid such a problem, the stellar neg- 
ativ e feedback is routinely invoked (see for r ecent works 
e.g.. lOkamoto et al.|[2Q05t iZavala et al.ll2008l though al- 
ternative physical processes have been also considered, 
e.g., Sommer-Larsen & Dolgov 2001; Governato et al. 
2004). In the same vein, it is well known now that the 
mass fraction of baryons in galaxies with respect to their 
halo masses is much smaller than the universal baryon 
fraction ff, = fi&/fio- The 'missing baryons' are com- 
monly assumed to have been expelled from galaxies by 
feedback effe cts, especially in the past, when the halos 
were smaller (|Dekel fe Silk 111986ft . 

In the absence of a general theory of SF and feed- 
back, several prescriptions have been proposed in order 
to model the corresponding subgrid physics in numeri- 
cal simulations at parsec scales or smaller, with the aim 
to properly extend SF/feedback influence into the larger 
(tens/hundreds of parsecs) resolved scales. Unavoidably, 
subgrid SF and feedback prescriptions in simulations in- 
troduce several free parameters. 

In simulations, some of the common criteria used for 
declaring a gas cell or particle as star forming are con- 
verging cold gas flows, a local Jeans instabili ty and/or a 
local gas density thre shold (tisf) (see e.g., iKatz et al~l 
119961: IKravtsov et~aTl |2003f ). The rate at which the 
gaseous mass element obeying the above criteria is con- 
verted into 'stars' is calculated using recipes which 
introduce some fr ee parameters like the SF efficieny 
constant 0. (e.g . . iKatzl 119921: INavarro fc White! 119931: 
IGerritsen fc Icke 1119971: ISpringel fc Hernquist II2003I ). In 
order to take into account the stochastic nature of gas 
conversion into stars, and also the fact that the mini- 
mum resolution element overlooks local variations in gas 
density, a probabilistic function around n^F is typically 
introd uced instead of a deterministic density criterion 
(e.g.. lStinson et al.ll2006l ). 

For the modeling of the feedback processes associated 
to the non-resolved regions, in which they start to oper- 
ate, two schemes are commonly used: individual gas ele- 
ments surrounding a stellar par ticle are given a velocity 
kick (kinema t ic fee d back: e.g.. [Navarro fc White! 119931 : 
lAbadi et all l200l iDubois fc Tevssier I l200l ~and/or 
are injected with t hermal energy and metals (e.g., 
IKravtsov et al.l [2003T ) . at the same time that, in some 
schemes, the radiative cooling is temporarly turn 
off to prevent the thermal energy from being ra- 
diated away before it is con verted into kinet i c en- 
ergy (adiabatic feedback; e. g .. IGerritsen fc Icke 1 119971 : 
iThacker fc Couchman I |200Ct iStinson et al.l 120061 ). Al- 
ternatively, in other schemes, the interaction between 
neighbor partic les belonging to different phases is inhib- 
ited by hand (Okamot o et al.l 120031 : iScannapieco et all 
2006). A third scheme, often applied, imposes an ef- 
fective gas equation of state in or der to describe the 
unresolved multiphase ISM (e.g., lYepes et al. I 119971 : 
ISpringel fc Hernquist I [2003). Hybrid schemes were also 
used: semi- analytical models of winds and starburts de- 
pendent on the gas physical conditions are included in the 



simulations (|Springel fc Hernquist 1 12003) . In all these 
feedback schemes there are free parameters, such as the 
amounts of momentum and thermal energy deposited 
into the gas elements, the time the cooling is switched 
off, t Q g, etc. 

The values for the free parameters in the SF and feed- 
back prescriptions are chosen in different ways. They can 
be defined a priori, based on constraints typically related 
to observations of the Galaxy ISM, and/or a posteriori, 
choosing the parameter values that make the evolution 
and properties of the simulated galaxy close to the ob- 
served ones. The prescriptions and parameters used also 
depend on whether the code is of the smoothed particle 
hydrodynamics (SPH) or Eulerian adaptive mesh refine- 
ment (AMR) kind. 

Because of the large diversity of models, strategies, 
and relevant parameters of SF/feedback implementa- 
tions, a comparison between them is not trivial. Al- 
though some implementations in high-resolution simu- 
lations have reached an amount of success forming rota- 
tional supported disk galaxies (see, fo r recent reviews, 
IMaver et all 120081 : IGibson et all 120091 and references 
therein), the task is far from being completed. It is not 
clear what are the quantitative differences between each 
subgrid model. Thus, a characterization of the various 
subgrid schemes for the ISM energy injection and dis- 
sipation is needed to gain some insight into their abil- 
ity/inability to form disk galaxies. 

In this paper, we explore the effects of the SF/feedback 
parameters on the structure and evolution of a low- 
mass galaxy in cosmological simulations using the Hy- 
drodynamics Adaptive Re finement Tree (ART) code 
(jKravtsov et al.lll997ll2003l ). We reach a maximum res- 
olution of about one hundred parsecs at redshift z = 1 
in order to be able to resolve the clouds associated with 
the SF; most of our analysis is made at this epoch and 
at z = 0.43, the last epoch reached by all the runs. Our 
nominal resolution extrapolated to z = is of 218 pc, 
which is among the highest resolutions for a cosmological 
galaxy format ion experiment usi ng the AMR/Eulerian 
technique (c.f. IGibson et al.l 120091) . 

We will focus our analysis on general evolutionary 
and physical properties/processes of the simulated galax- 
ies, namely: (i) the evolution of structural and dy- 
namical properties (e.g., galaxy stellar structure, cir- 
cular velocity curve decomposition), (ii) the struc- 
ture/thermodynamics of the ISM (density and temper- 
ature probability distribution functions, vertical filling 
factors of the different phases of the ISM, gas fractions), 
(iii) the dependence of the global SF history and the 
SFR on gas density, and (iv) global quantities like the 
baryon mass and angular momentum fractions. Our aim 
here is not to simulate galaxies that can be compared to 
observations, but rather to determine the effects of the 
parameters of the SF/Feedback prescriptions. Moreover, 
we will use our results to point out some potential diffi- 
culties of simulated low-mass galaxies in the context of 
the ACDM cosmogony. 

The plan of the paper is as follows. The code and the 
SF and feedback prescriptions used for the simulations 
are described in §2. In this section, the cosmological sim- 
ulation and the different runs varying the SF/feedback 
parameters are also presented. The results of the various 
runs are given and compared with each other in §3. In 
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§4 we present a summary of the main results obtained 
in the paper accompanied by interpretations and discus- 
sions. Some potential vexing problems of any ACDM- 
based model of low-mass galaxy formation and evolution 
arc hihglightcd in §4.3. Finally, our conclusions are given 
in §5. 

2. NUMERICAL EXPERIMENTS 
2.1. The Code 

The numerical simulations described here were per- 
formed using the hydrody namics + N-body Adaptiv e 
Refinement Tree code ART (|Kravtsov et al.lll997l [2003h . 
Among the physical processes included in the code are 
the cooling of the gas and its subsequent conversion into 
stars, stellar feedback, self-consistent advection of met- 
als, an UV heating background source, etc. 

The cooling and heating rates incorporate Comp- 
ton heating/cooling, atomic and molecular cooling, 
UV heating from a cosm ological background radiation 
(jHaardt fe Madau Ill996h . and are tabulated for a tem- 
perature range of 10 2 < T < 10 9 K and a grid of den- 
sities, inetallj^hks, and_redshifts using the CLOUDY 
code (jFerland et al. II1998L version 96b4). Gas is con- 
verted into stellar particles according to a prescription 
described and discussed below. Stellar particles eject 
metals and thermal energy through stellar winds and 
type II and la supernovae explosions. Each supernova 
dumps 2 x 10 51 ergs of ther mal energy and ej e cts 1 .3M^ 
of metals. For the assumed iMiller fc Scalo I ()1979l ) ini- 
tial mass function, IMF, a stellar particle of 10 5 M Q 
produces 749 type II supernovae. For a more detailed 
discussion of the processes i mplemented in the code, see 
iKravtsov etHI ([20031 [2001 . 

2.2. Star formation and feedback 

In ART, SF is modeled as taking place in the coldest 
and densest collapsed regions, defined by T < T$f and 
Pg > Psf, where T and p g are the temperature and den- 
sity of the gas, respectively, and psf is a density thresh- 
old. A stellar particle of mass m* is placed in a grid 
cell where these conditions are simultaneously satisfied, 
and this mass is removed from the gas mass in the cell. 
The particle subsequently follows N-body dynamics. No 
other criteria are imposed. 

Although nsF, the hydrogen number density threshold 
corresponding to psf, is in principle a free parameter of 
the simulations, one can make an "educated guess" as to 
which values should be considered most realistic, by not- 
ing that SF occurs almost exclusively in giant molecular 
clouds (GMCs). Thus, it is natural to require that the 
density in a star-forming cell should be at least such that 
the cell's column density equals the threshold value for 
molecular Hydrogen and CO formation, 7V ^ 10 2 1 cm~ 2 
(jFranco fc Cox|[l98l Ivan Dishoeck fc Black) [l9l8h . For 
a cell of 150 pc, this translates into a number density 
Usf ~ 5cm -3 . In our models, tt-sf is varied around this 
value from 0.1 to 50 cm -3 . We keep Tsf fixed at 9000 
K. 

The stellar particle mass, m„, is determined through 
a simple two-step process. A first estimate of this mass, 
labeled mi , is made by assuming that m* is proportional 
to a power law of the gas density, . In this work we use 
a = 1.0. Thus, in this step, collisionless stellar particles 



can be created with dto, where C* 

is a constant that measures the "efficiency" with which 
gas is locally converted into stars during dto, the time 
step of the coarsest grid, which is also the time interval 
between successive checks by the code for the formation 
of new stellar particles; m gas is the gas mass in the cell. 
For low values of C*, for example C* = 2.5 x 10~ 10 yr" 1 
(the fiducial value in the code), the formula above pro- 
duces small mi values. However, this makes the num- 
ber of stellar particles in the simulation practically in- 
tractable. For this reason, a stellar particle mass limit is 
introduced in the SF set of parameters, labeled m*,u m . 
The SF algorithm then chooses the maximum between 
mi and m*.i; m , which we call mi- In the second step, the 
algorithm takes the minimum between mi and 2/3 of the 
gas mass contained in the cell and assigns this mass to 
the stellar particle, in order to prevent total exhaustion 
of the gas contents of the cell. Since the stellar particle 
masses are much more massive than the mass of a star, 
typically 10 4 - 10 5 M , once formed, each stellar particle 
is considered as a single stellar population, within which 
the individual stellar masses are distributed according to 
the Miller & Scalo IMF. 

The simulations reported here are of two kinds: those 
with a "deterministic" SF prescription (the "D" runs), in 
which stellar particles are always created once the den- 
sity and temperature conditions are satisfied, and those 
with a "stochastic" or "random" SF prescription (the 
"R" runs ), in which stellar particles are created in a cell 
with a probability given by 

!0 if P g < Psf; 
1 if/3p S F<Pg- 

where j3 is a free parameter, to which we refer as the 
probability-saturation threshold. Thus, regions with 
higher densities have higher probabilities to host SF 
events. This stochastic prescription allows for the possi- 
bility of forming stars in regions of low average density 
(assuming we have chosen a low tisf parameter) , in which 
isolated dense clouds would not be resolved. 

The runs can also be divided in two classes according 
to the stellar feedback implementation. In the first class 
of runs, all the energy from stellar feedback is dumped 
onto the star-forming cel l in the form o f heat , as de- 
scribed, for example, in IKravtsov et al.l (|2003l ). This 
thermal feedback suffers from the well-known overcool- 
ing problem in low-resolution simulations. This problem 
can be traced back to ins ufficient numerical resolution 
(|Ceverino fc Klypin II2009T) , and happens when the char- 
acteristic expansion time in a cell heated by a stellar par- 
ticle (given by the sound crossing time across the cell) is 
much longer than the cooling time at high densities and 
low resolutions (large cell size). In this case, in a single 
timcstep, the gas cools before it has time to expand. The 
resolution condition to avoid this problem would then be 
to require that the sound crossing time across the cell be 
shorter than the local cooling time. This, however, can 
require impossibly high resolutions with presently avail- 
able computational resources. 

A fr equent strategy to o v ercome this problem 
(e.g., iGerritsen fc Icke I 119971: iThacker fc Couchman I 
120001: iSommer-Larsen et al. I 120031: iKeres et al. I 12005 
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iGovernato et al.ll2007f) is to turn off the cooling in a cell 
for some time after a stellar particle is produced there, in 
order to allow for the gas to expand before it cools. Al- 
though this procedure has been criti cized as being some- 
what arbitrary and unrealistic (e.g., iCeverino fc Klypin I 
2009), we consider that it is physically motivated, as it is 
even more unrealistic to prevent the expansion of heated 
regions due to the inability to resolve the scales that 
would allow the physically correct behavior of these re- 
gions. In the present paper we therefore adopt this strat- 
egy, and test its effects on the evolution of the forming 
galaxy. Specifically, we turn off the cooling in a cell for 
40 Myr after a SF event occurs there in some of the runs 
(see Table [J). 

2.3. The numerical method 

We present here several simulations of the evolution 
of a galaxy with a total (baryon + dark matter) mass of 
about 7 x 10 10 h~ 1 M Q at the present day. The analysis is 
made at two epochs, z — 1 and z — 0.43. The latter cor- 
responds to a time just before the last major merger. The 
simulations are run in a ACDM cosmolgy with fio = 0.3, 
Qa = 0.7, and tth = 0.045. The CDM p ower spectrum is 
taken from lKlypin fc Holtzmanl (|1997[ ) and it is normal- 
ized to erg = 0.8, where os is the rms amplitude of mass 
fluctuations in 8 ft, _1 Mpc spheres. 

To maximize the resolution efficiency, we first perform 
a low-resolution N-body simulation with 128 3 dark mat- 
ter (DM) particles in a periodic box of 10/i _1 Mpc on a 
side. At the start of the simulation, the box is initially 
covered by a mesh of 128 3 cells (zeroth level cells). A 
spherical region, centered on a randomly selected halo of 
mass 2.6 x 10 11 /i~ 1 M Q and with a radius equal to three 
times the virial radius (i? v ir) of the halo, was then chosen 
at z — and its corresponding Lagrangian region, identi- 
fied at z = 50, was re-sampl ed with additional small-scale 
modes (|Klypin et al. ||200TI ). The virial radius is defined 
as the radius that encloses a mean density equal to A v i r 
times the mean density of the universe, where A v ; r is a 
quantity that depends on D,q, Q\ and z. For example, for 
our cosmology A v j r (z = 1) = 203. The number of DM 
particles in the high-resolution zone is about one million 
and the mass per particle (mdm) is 5.3 x 10 5 /i~ 1 M Q . 

In ART, the initially uniform grid is refined recursively 
as the matter distribution evolves. The criterion chosen 
for refinement is based on gas or dark matter densities. 
The cell is refined when the mass in DM particles exceeds 
1.3(1 — fb)m p or the mass in gas is higher than 13.0/bm p , 
where m p = md m / (1 — /&)■ For the simulations presented 
in this paper, using multiple dark matter particle masses, 
the grid is always unconditionally refined to the third 
level, corresponding to an effective grid size of 512 3 . On 
the other hand, the maximum allowed refinement level 
'max was set to 9. At z = 1, the number of grid cells, at 
all levels of resolution, is about 11 million. In particular, 
at l max the number of cells is ~ 5 x 10 5 , corresponding 
to a resolution of 218 comoving parsecs. 

In this paper, we focus our study on the third most 
massive halo appearing in the re-sampled simulations at 
z = 1 because it turns out that it evolves from z ~ 2 
to z ~ 0.4 in an nearly isolated fashion. This allow us 
to observe the effects that the SF and feedback prescrip- 
tions have on the structure and evolution of the embed- 
ded galaxy without having to deal with complicated dy- 



namical interactions. This halo contains within its virial 
radius ~ 10 5 DM particles at z = 0. In the remain- 
der of the paper, we study the properties of this system 
as a function of the SF/feedback prescriptions and the 
corresponding parameter values. 

2.4. The runs 
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Fig. 1. — Stellar particle mass distributions measured at z = 1.0 
for all simulations presented in Table [T] The distribution depends 
on the values of the SF parameters m, i; m , C*, and nsp. 

In Table [T] we list the parameters described above for 
the various runs we performed. The first column gives 
the name of the run, denoted by a capital letter followed 
by a number, where the former indicates whether the 
SF model belongs to the D or R sequence (see §2.2| . In 
columns (2), (3), and (4) we show the values of m^um, 
nsF, and C*, respectively. Column (5) indicates whether 
cooling is always on or it is turned off for 40 Myr after 
SF events in the simulation. The parameter /?, defined 
above (see eq. [IJ, is shown in column (6). We also ran 
an adiabatic simulation. This is denoted in Table Q] by 
the name Ad (see last row). 

The SF scheme and the different values of the SF pa- 
rameters TO*,ii m , nsF, and C* produce different distribu- 
tions for the stellar particle masses, mJ3- Figure [T] shows 
this distribution at z = 1.0 for each of the runs presented 
in Table [TJ Some features can be highlighted from this 
figure: (a) the distributions of the R sequence of runs 
are wider than the corresponding distributions of the D 
sequence. This just reflects the fact that the recipe for 
SF for the former case is probabilistic; that is, stellar 
particles can be produced from a wide range of cold gas 
density values, (b) In most runs the distributions have 

1 These masses refer to the initial masses with which stellar par- 
ticles were born. At later times, the masses are always lower than 
the initial (original) masses because stellar particles loose mass due 
to stellar winds and supernova events 
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TABLE 1 

Parameter simulations 



Model m„,iim a n SF b C* c Cooling d /3 C 

(10 4 'M Q ) (cm- 3 ) (2.5 X 1CT 10 yr" 1 ) (10 6 yr) 



Dl 1.0 50.0 1.0 on off 

D2 8.0 6.25 1.0 on off 

D3 8.0 6.25 1.0 40 off 

D4 oo 6.25 none 40 off 

Rl 1.0 1.0 667.0 on 100.0 

R2 1.0 0.1 500.0 40 20.0 

R3 1.0 1.0 500.0 40 20.0 

R4 1.0 6.25 500.0 40 20.0 

Ad off off off off off 



Note. — Runs are denoted with the capital letters D or R followed by a 
number. D stands for deterministic and R for random. The adiabatic run is 
denoted by Ad. This simulation does not include radiative cooling. 
a Stellar particle mass limit. 
b Threshold gas density for star formation. 
c Star formation efhcicny constant. 

d Cooling is always on or turned off for 40 Myr after SF events. 
e Probability-saturation threshold. 



a spike at their corresponding rn*fi m value, (c) From se- 
ries D, run D4 is the one with the largest values of m*. 
In this run, where TO*j; m is set to infinite, m* is always 
equal to 2/3 the gas mass in the cell, (d) Run R2 has 
the widest distribution and, along with run R4, it shows 
the largest values of m*. Below, we briefly discuss the 
characteristics of each run. 

Run Dl has m*fi m — 10 4 M©, ngp = 50, and the 
cooling is always on (see Table [1]). The masses of stellar 
particles vary from about 10 4 M© to 6 x 10 4 Mq and peak 
at 10 4 M©. This run has the highest density threshold 
and, along with run R4, the most peaked circular velocity 
curve (see £ )3.2[) . 

Run D2 was performed to test the effect of a lower 
density threshold, although this necessarily implies that 
more stellar particles will form. We thus increased the 
value of m* l iii n from 10 4 M Q to 8 x 10 4 M© in order 
to keep the number of stellar particles computationally 
tractable. Both of these models have a low C* value 
but m»j; m is higher in model D2. This means that, in 
practice, in this model, in the first step of the SF pre- 
scription, 77i2 = ™*,iim always (see £|2.2j) . Note, however, 
that if mi — m^um then the final chosen mass of the 
stellar particle cannot be higher than m»ji m . Thus, in 
model D2 the masses of stellar particles are restricted to 
be lower than 8 x 10 4 M©. 

Run D3 is similar to D2 but here we turn off the cooling 
for 40 Myr after SF events. Run D4 is similar to D3 but 
here m».ii m = oo, which means that cells always form 
stellar particles with a 2/3 efficiency factor (see i)2.2[) ; 
that is, in any SF event, 2/3 of the gas mass in the cell 
is converted into stars. The masses of the stellar parti- 
cles vary from about 10 4 M© to 4 x 10 5 M©. This model 
is expected to have a higher impact, as compared with 
models D2 and D3, on the ISM properties of the simu- 
lated galaxy due to the formation of more massive stellar 
particles, which are concentrated sources of high energy 
injection. 

Aside from m*ji m and nsF> the R sequence of models 
include the j3 parameter (see eq. [T]). In model Rl this 
parameter is 100 while in the rest it is 20; model Rl also 
differs from the rest in that it includes "runaway stars" 



(|Ceverino fc Klvpin~ll2OO90 . stellar particles able to inject 
energy to the ISM in regions of low-density gas. In this 
simulation, which uses a different version of the code, a 
random velocity, drawn from an exponential distribution 
with a characteristic scale of 35 km s _1 , is added to the 
inherited gas velocity to 30% of the stellar particles (see 
Ceverino & Klypin 2008 for more details). We used a 
characteristic velocity twice higher than the one used by 
Ceverino & Klypin, whose value was motivated by obser- 
vations, to take into account the fact that our simulations 
are twice less resolved. The masses of stellar particles in 
the R runs vary from about 10 3 M© to 10 6 M©. Runs 
R2, R3 and R4 differ only in the value of tisf- All have 
«T-*,iim = 10 4 M©, f3 — 20.0, and the cooling switched off 
for 40 Myr after SF events. They were run to investigate 
the effect of varying tisf on the structural and evolution- 
ary properties of the simulated galaxy in the stochastic 
SF scheme. 



3. RESULTS 

As mentioned above, our study focuses on the third 
most massive halo found at z = 1 because this halo 
evolves without significant mergers, at least in the red- 
shfit range from z ~ 2 to ~ 0.4. Most of the analysis of 
the embedded galaxy is performed at two epochs: z = 1 
and 0.43; for the cosmological model used here, the time 
interval between these two redshifts is 3.2 Gyr. 

Some global properties of the halo in the adiabatic sim- 
ulation, at different redshifts, are presented in Table [5J 
Column (4) shows the concentration of the halo as de- 
fined by c v i r — i? V ir/?~J3- The virial mass of the halo and 
the specific angular momentum are given in columns (2) 
and (6), respectively. The baryonic fraction within i? v j r 
is presented in column (3) (see also Table . Column 
4) shows the spin p arameter of the halo as defined by 
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Fig. 2. — Spatial distribution of stellar particles for runs D4 and R2 at z = 1 and z = 0.43; the disks, when present, are shown edge 
on. Scale is comoving h _1 kpc (solid line measures 10 ft - x kpc). Particles are color coded by age with a rainbow palette with the bluest 
particles being the youngest ones. At z = 1 in run D4 a young stellar disk can be clearly distinguished while in run R2 the young stellar 
particles are distributed in a spheroid. 



where J is the angular momentum inside the halo virial 
radius, and V V i r is the circular velocity at i? v ir- More- 
over, column (2) of Table [3] reports the virial masses of 
the DM halo in the other runs. At z = 1, M V - 1T varies be- 
tween 4.7 and 5.0 x 10 10 M Q , while at z = 0.43, it varies 
between 5.4 and 6.6 x 10 10 M©. In the following, we will 
discuss the dependence of the galaxy morphology, stellar 
populations, dynamics, ISM properties, and SF history 
of the studied galaxy on the relevant parameters of the 



SF and feedback prescriptions. 

3.1. Morphology and structure 

The different combinations of the parameters used in 
the SF and feedback prescriptions, described in §2.4, pro- 
duce a large diversity of galaxy morphologies and stellar 
populations. As an example, in Fig. [5] we show the stel- 
lar spatial distribution for runs D4 and R2 at z = 1 
and z = 0.43, with the disk, when present, seen edge-on. 
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TABLE 2 

Adiabatic Halo Global Parameters 



M vir / 9a ((Rvir) a <W A J/M b 

(10 10 M Q ) (km s^ 1 kpc) 



0.0 8.6 0.14 18.4 0.044 285.1 

0.4 5.7 0.14 13.4 0.049 212.4 

1.0 4.8 0.14 8.0 0.060 202.8 



a Baryon fraction within the virial radius of the halo. 
b Specific angular momentum inside -R v i r . 

Stellar particles are color coded according to their age in 
a rainbow scale, with the reddest particles being oldest. 
Ages range between about 6 Myr to 5.5 Gyr (z = 1.0) and 
12 Myr to 8.7 Gyr (z = 0.43). Unlike run D4, in which 
a young stellar disk embedded in an older thick disk can 
be appreciated, the young stellar population of run R2 at 
z = 1.0 is distributed in a "puffy" spheroidal structure. 
At z = 0.43, 3.2 Gyr since z = 1, the once young stellar 
disk of run D4 has aged, while in run R2 a young stel- 
lar disk has formed. This illustrates the large effects the 
SF/feedback prescription has on the s tructure of the an- 
alyze d galaxies (see also, for example. iScannapieco et al.l 
l2008h . 

In column (5) of Table [31 we report, for z = 1.0 and 
0.43, the galaxy stellar effective radius, R e , defined as 
the radius where half of the central galaxy stellar mass 
M* (disk and spheroidal components) is contained. We 
find that the galaxy stellar structure is systematically 
larger as usf decreases (sequences Dl — > D2 for the de- 
terministic model, and R4— > R3— > R2 for the random 
model). This radius is particularly small for those runs 
where feedback is inefficient, Dl and D2; when the radia- 
tive cooling is artificially switched off allowing for ther- 
mal pressure-driven feedback, the radius increases by a 
factor of ~ 2 (from run D2 to D3). However, R e seems 
to be insensitive to the stellar particle mass distribution 
(compare runs D3 and D4). 

Figure |3j show the stellar surface density profiles along 
the disk plane, averaged azimutally in concentric cylin- 
drical rings, for all the runs at z = 0.43; radii are nor- 
malized to their corresponding R e . In most cases, the 
surface density profiles can be described approximately 
by two exponential laws: an inner nearly exponential 
profile associated to a disk+bulge system and a more ex- 
tended outer exponential profile associated to an older 
spheroidal structure (stellar halo). It is interesting that 
the outer profiles of the R models tend to be higher and 
more extended than those of the D models. This is prob- 
ably because models with a stochastic SF scheme favors 
star formation in low gas density environments, helping 
the formation of an exponential stellar halo. 

The measured effective stellar surface densities of the 
various runs, T,m = A/*/27ri?g, lie in the higher density 
side but within the range of observational determinations 
at z ss 0.4 and z ss 1; in other words, the simulated galax- 
ies have smaller R e for their masses t han the average rad ii 
inferred from observations (see e.g., [Bardcn ct al. 2005). 

In order to quantify the relative contribution of the 
spheroidal and disk components to the total galaxy stel- 
lar mass we compute the ratio, e, between the z com- 
ponent of the angular momentum of the stellar particle 
and the angular momentum expected for a circular orbit, 




\ 



t i I i I i I i I i x . r , I , i , i , I , I j 
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Fig. 3. — Stellar surface density profiles for the D (left panel) 
and R sequence of runs at z = 0.43. 

j c = RV C (R). We roughly define the mass of the spheroid 
as twice the mass of all stellar particles with e < 0. The 
mass of the stellar disk is then estimated as the differ- 
ence between the stellar mass reported in Table [3] and 
the mass of the spheroid. In all runs, the spheroidal stel- 
lar component is significant (between ~ 30% and ~ 60% 
of M* at z = 0.43) and with old/intermediate ages. The 
runs with the best defined stellar disks are Rl, D3, and 
R3. The spheroid-to-total stellar mass ratio for these 
runs is ~ 30 — 35%; the complement is roughly the disk 
component. It should be mentioned that run Rl fol- 
lows a m odel constructed to reprod uce Milky- Way-sized 
galaxies (jCeverino fc Klypin~| [2009() . 
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Fig. 4. — Circular velocity profile decompostions for the D se- 
quence of runs at z = 1: gas (short-dashed line), dark matter 
(dotted line), stars (long-dashed line), and total (solid line). Only 
the total profiles are shown also for the epoch z = 0.43 (thick dot- 
dashed line). Runs in which cooling is always on, Dl and D2, show 
very peaked profiles, being run Dl the most extreme. On the other 
hand, run D3 and specially run D4, show relatively flat profiles. 



3.2. Circular velocity decomposition 

Figures H] and [S] show the total circular velocity profiles 
for the R and D models described in Table Q] at z = 1.0 
(solid line) and 0.43 (dot-dashed line). The circular ve- 
locity is defined as V c = ^JGM (R) /R, where M(R) is 
the total mass, or the mass of a certain galaxy compo- 
nent (dark matter, gas or stars), contained within radius 




Fig. 5. — As in Figure[4]but for the R sequence of runs. Here we 
have also added the rotational velocity of the cold gas (T < 10 4 K) 
that is inside the disk and plotted with a solid line with symbols 
(squares). Runs R2, R3, and R4 differ only in the value of ngp and 
it increases as we go from R2 to R4. A clear trend of increasing 
V c and steepening of the profile is seen. Runs R2 and Rl have a 
relatively fiat circular velocity profile. 



R. The circular velocities due to the various mass com- 
ponents are plotted only for epoch z = 1.0. 

In both, the deterministic (D2 — > Dl) and the stochas- 
tic (random) runs (R2— >• R3— !• R4), we see that an in- 
crease in nsF translates into more centrally concentrated 
stellar and gas structures, which produces a more peaked 
total circular velocity profile. The halo also becomes 
more concentrated, but this is a consequence of the bary- 
onic gravitational drag. As expected, for our spatial 
resolution, the contribution of the gas component to V c 
increases as hsf increases. This trend might change if 
tens of parsec resolution is reached. The lack of efficient 
feedback produces stellar disks that are too concentrated 
(runs Dl and D2); when local cooling is switched off tem- 
porarily, keeping all other parameters fixed, V C (R) flat- 
tens significantly (run D2 to D3). 

It is interesting to note that run R2 has a mild de- 
clining circular velocity curve. However, this run has no 
disk at z = 1, as seen in Figure [2j Nevertheless, a disk 
forms later on and its V c becomes sliglthly more peaked. 
Runs D3 and D4 also have mild circular velocity curves 
at z = 1, but, unlike run R2, these runs do have well 
defined young stellar disks at this redshift. As time goes 
on, the circular velocities of these models become more 
peaked and the galaxies more concentrated (see Figure 
3]). Interestingly enough, such an evolution is much less 
pronounced for runs with the random SF prescription 
(see Figure [5]). 

Lastly, it is important to say that the average rota- 
tional velocity in the mid-plane of the disk docs not trace 
the galaxy circular velocity in our simulated galaxies (see, 
for example, Figure [5]). In some of our models, the tan- 
gential velocity profile slowly increase s with radius as it 
is observed in low-mass galaxies (e.g.. lPersic et al.lll996l 
but see Swaters et al. 2009). Notice, however, that 
at the current resolution our galaxies do not present a 
central bar. Similar situations in which rotation curves 
do not trace the circular velocity have been discussed 



bv lValenzuela et al. I (|2007t ) and iRhee et al. I (|2004f >. In 

these cases the difference is attributed to pressure gradi- 
ents and non-circular motions. 

3.3. Star formation history 

Figure [5] shows the SF history (SFH) for the various 
runs, defined as the instantaneous SFR as a function of 
time, and computed for each run using the last snapshot 
recorded. Specifically, in each data dump, in addition 
to the positions and velocities for all stellar particles, we 
also save the time at which they formed, their masses 
(initial and present), and metallicities due to supernovae 
of types la and II. We add up the (initial) masses of the 
stellar particles formed during a certain time interval, 
which we take as 0.1 Gyr, and divide it by this time, 
to obtain the "instantaneous" SFR. Since the amount 
of stellar mass outside the galaxy is small during the 
epochs analyzed (< 3%), the SFH computed inside i? v j r 
is basically that of the galaxy. The simulations were run 
until z = 0.43, whichs corresponds for the cosmology 
used here to a cosmic age of 8.9 Gyr. 

In all runs, the SFHs are characterized by a fast in- 
crease lasting in most cases 1.5-2.0 Gyr followed by a 
nearly exponential decreasing phase. Runs with smaller 
values of usf show a more gradual decrease in general 
and the peak of the SFH is shifted to later epochs (se- 
quences R4— > R3— > R2, and Dl— > D2). In the case of 
the R sequence, the increasing phase is also more grad- 
ual and the amplitude of the maximum decreases with 
decreasing usf- Run R2, which has nsF = 0.1 cm -3 , 
exhibits the most gradual SFH initial increase and the 
latest SFH maximum (t ~ 3.5-4 Gyr). Runs R3 and R4 
are those with the fastest SFR decline with time. 

The SFH is also stronlgy affected by the stellar ability 
to heat the surrounding gas. The amplitude of the initial 
maximum in run D3, in which the cooling is temporarily 
switched off, is lower, and its subsequent decrease is more 
gradual, as compared with run D2. The SFH of the 
former run is much burstier than the one of the latter run, 
with periods of very low SF activity followed by intense 
bursts. This bursty behaviour is even more pronounced 
in run D4, a run similar to D3 but with more massive 
stellar particles (compare also the SFHs of runs Dl and 
D2). Notice also that in general runs with the random SF 
prescription tend to have a less bursty SFH than those 
with the deterministic one. 

3.3.1. Relationship between star formation rate and gas 

density 

Figure [7] shows the average SF rate (SFR) surface den- 
sity, Ssfr, versus the gas surface density, H g , for all the 
runs at z = 1 and z = 0.43. The surface SFRs are com- 
puted using the surface densities of stellar particles with 
ages smaller than 4.0 x 10 7 yr. For the gas surface density, 
only gas with T g < 10 4 K is considered. Cylindrically- 
averaged values of both quantitites with a constant radial 
bin of 0.5i? e wide and 0.5 /i _1 kpc proper height are ob- 
tained for each model and each epoch. For each run, a 
number of points are plotted with this number depend- 
ing on the extent of the simulated galaxy. The symbols 
corresponding to each run are given in the figure caption. 
The symbols connected by the dotted and solid lines cor- 
respond to z = 1.0 and z — 0.43, respectively. The large 
symbols correspond to averages of Esfr and S g inside 
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Fig. 6. — SFR histories for all runs presented in this paper. Each 
panel shows two runs labeled inside panels. In the X axis runs the 
cosmic time. The last time shown corresponds to z = 0.43. For 
given run, the SFR as a function of time is computed using the 
snapshot at z = 0.43 by tracing back the recorded formation times 
and initial masses of the stellar particles that are inside the virial 
radius. 

3i? e . Notice that surface gas densities move to the right 
(left) if we take a longer (shorter) cylinder. 

As expected, as we go to smaller nsF values (sequences 
R4— > R3— > R2, and Dl— > D2), the points move down and 
to the left side of the £sFR~£ ff diagram. Stellar feed- 
back (see for example, from run D2 to run D3), does not 
impact significantly the EsFR-Eg diagram; at z — 0.43 
these runs occupy a similar locus in this diagram. As 
time goes on, while E g decreases slightly, Ssfr decreases 
strongly for most of the cases. In general, for a given sim- 
ulation, the correlation of Esfr with E g , both along the 
radius and at two epochs, is steep and becomes steeper 
as S g decreases. 

3.4. Properties of the interstellar medium 

3.4.1. Probability density functions (PDFs) 

Figures [5] and M show the normalized probability den- 
sity function, PDF, of gas number density for the D and 
R series, respectively. Thin solid and dotted lines refer to 
measurements at z = 1. The latter show the PDFs com- 
puted using all cells in the maximum level of refinement, 
I = inm, while the former represent the PDFs computed 
using the cells contained within a cylinder with a radius 
of 3R e and height of 0.6 ft, _1 kpc proper, centered on and 
oriented with the disk. The normalization volume, Vt t, 
differs in both cases, being smaller for the cylinder. The 
thick dot-dashed line refers to the PDF measured for the 
gas within the cylinder at z = 0.43. The plane of the disk 
is defined as the plane normal to the angular momentum 
vector of the cold gas located in an inner cylindrical shell 
of radii of 0.5 and 3 /i _1 kpc comoving. We have varied 
the inner and outer radii of the cylinder and found negli- 
gible variations in the measured quantities. Notice that 
PDFs computed using the cells at I = Z max include a sig- 
nificant fraction of cells that contain low-density warm 
and hot gas. This is due to the fact that these low-density 
gas cells are kept refined because the dark matter density 
cell refinement criterium is satisfied. 




1 10 100 1000 



Fig. 7. — Radial SFR surface densities, Egpri, versus radial gas 
surface densities, S ff , for each one of the runs at two epochs. The 
points of a given run are plotted with the same symbol and con- 
nected with (red) dotted lines for the epoch z = 1.0 and with solid 
(black) lines for z = 0.43. The large symbols show the correspond- 
ing average value of Ssfr and S 9 of each run at the two epochs 
(see the text for more details on how the radial and average Sgpj^ 
and S 9 were calculated ). The symbol code for each run is as 
follows: Dl (open square). D2 (skeletal, four-end), D3 (open star, 
four-end), D4 (solid square), Rl (open pentagon), R2 (skeletal, 
five-end), R3 (open star, five-end), a nd R4 (solid penta gon). The 
solid line is the fit to observations by Kcnnicutt] l|1998I V 



-3 



> -Z 



-3 



1 

- Dl 


1 i 1 i 1 i 1 i 1 i 

ft A 
r 1 V \ 

L,:--= J AH : k 




• D2 b 






j 


-, 


- 
_ 


- ] ' n 

Jr 

ni * 


I. 

i 

H r 


)! 








A 

_ D4 i \ 






/ 

r J / 1 
\ f 


, 

1 

1. 
1 

h 




/ / / i r 

J I J 
f J 


r 


1 


! 

r, i J i,i 


1 

ill 




+ I r 

f r i , i 




-4 


-3 -2 -1 
log n H (cm 3 ) 


1 £ 




-4 -3 -2 -1 
log n H (cm 3 ) 


1 2 



Fig. 8. — Probability distribution functions of gas density, hh, for 
the runs from series D. The rijj PDF's, the fraction of total volumen 
occupied by gas in a given density logarithmic bin, were computed 
at z = 1 or 0.43 using all cells in the maximum level of refinement 
(dotted lines) or using the cells contained in a cylinder of radius 
3-Re and height 0.6 h~ kpc proper (solid and thick dot-dashed 
lines) centered on the disk. At z = 0.43 only the PDFs within the 
cylinder (thick dot-dashed lines) are shown in the figure to avoid 
line crowding. The PDF's computed using the maximum level of 
refinement cover a bigger volume including a greater normalized 
fraction of low-density gas. 

Figures ITTJI and ITTI show the corresponding PDFs of the 
gas temperature; the line code is as in Figures[5]and[S] As 
was already noticed, the volume occupied by the cylinder 
is smaller than that of the cells at the maximum level of 
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Fig. 9. — As in Figure l8l but for the models of the series R. 

refinement. In particular, the temperature PDF of run 
Dl computed over the cylinder at z — 1 contains only 
cells with T < 10 4 K (see upper- left panel of Figure [TU]) : 
that is, the PDF of this model is only sampling the well- 
defined cold gaseous disk. Instead, the PDFs computed 
using the cells at I = Z max contain a significant fraction 
of warm-hot /hot gas. 

Figures IMTT1 show that at z = 1.0 (the situation is 
similar at z = 0.43, albeit more extreme) all models ex- 
cept R2 have very different PDF's when calculated over 
the cylinder around the mid-plane than when computed 
using the grid cells at I = l ma x- The main difference is 
the relative absence of hot, low-density gas in the mid- 
plane: this gas occupies a very small fraction of the disk 
volume. In particular, we see that in runs D3 and D4, 
the maximum temperatures in the disk are T ~ xq 5 5-6 
K. On the other hand, we see that the temperature PDF 
of the finest grid of run D4, in particular, contains a large 
amount of hot gas, and that its distribution reaches tem- 
peratures of up to 10 7 K. This suggests that the hot gas 
produced by stellar feedback in the disk of run D4 be- 
comes buoyant and leaves the disk. Indeed, run D4 even 
exhibits a wind (cf. ^3~4~2|) . 

Run R2 is the only one at z = 1 for which the density 
and temperature PDF's are very similar when measured 
for the mid-plane cylinder and for the finest grid. This 
can be understood because the galaxy in this run has 
a more spheroidal shape at this redshift. However, at 
z = 0.43 (thick dot-dashed line) a thick disk is already in 
place (see Figure^]) and the PDF's measured for the mid- 
plane cylinder and for the finest grid start to deviate from 
each other, with the low-density hot disk gas buoyantly 
escaping to the halo. 

Let us analize the PDFs calculated over the cylinder 
around the rotating mid-plane. We see that as nsF is 
made smaller (sequences R4— ^ R3— !• R2, and Dl— > D2), 
the density and temperature PDFs become broader. The 
main effect of the broadening is in the direction of in- 
creasing the amount of low-density and warm-hot/hot 
gas. In the extreme case of model R2, at z = 1 a second 
peak appears in the low-density and high-temperature 
sides of the corresponding PDFs (see Figures [9] and [TT]) . 



The maximum gas density reached by a run also depends 
on the value of n^F\ the lower tisf, the lower the maxi- 
mum Tin value. 

Run D2, in comparison to D3 fFigureslSlandfTU)). shows 
a higher fraction of low-density, hot gas in both the cylin- 
der or in the highest-refined grid cells (dotted lines). We 
see here that a stronger feedback (run D3) produces a 
less dominant hot gas phase, very likely because part of 
this hot gas is expelled from the mid-plane of the disk. In 
any case, the galaxy in run D2 is too compact (see Table 
[3]) and shows a rather thick young stellar disk, with its 
gaseous structure resembling that of run R2. The region 
in which we compute the PDFs contains a mixture of cold 
and warm-hot/hot gas, heated very likely by shocks. 

From run D4 to run D3 (Figures [8] and [TO]) , the PDFs 
of density and temperature both within the disk at z — 1 
(thin solid lines) and in the cells at the maximum refine- 
ment level (dotted lines) show a decrease in the lower 
density and higher temperature regimes. Therefore, by 
allowing for more massive stellar particles (a more effi- 
cient conversion of gas to stars), the galaxy produces a 
higher fraction of low-density and hot gas. 
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Fig. 10. — Probability distribution function of gas temperature, 
T, for the runs from series D. The PDF's of temperature, the frac- 
tion of total volume occupied by gas in a given temperature loga- 
ritmic bin, were computed as explained in the text and in Figure 
\E\ The same line code as in Figure [8] is used here. 

3.4.2. Vertical fillmg-factor profiles 

The spatial distribution of the gas in the runs can also 
be analysed through the filling factors of the gas at var- 
ious temperature ranges. The vertical profiles at z = 1 
are shown in Figures [12] and [13] We compute the filling 
factor over cylinders that have the same radius as the one 
we use to compute the PDFs; i.e., 3i? e of the correspond- 
ing simulated galaxy. We define four regimes according 
to the (azimuthally averaged) gas temperature: cold gas 
(T < 3000 K, solid lines), warm gas (3000 < T(K) < 10 4 , 
dotted lines), warm- hot gas (10 4 < T(K) < 10 5 5 , short- 
dashed lines ), and hot gas (T > 10 5 ' 5 K, long-dashed 
lines). 

Figures [12] and [13] show that as ngp decreases (se- 
quences R4— > R3— > R2, and Dl— > D2) and/or m*ji m 
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Fig. 11. — As in Figure llOl but for runs of series R. 



increases (sequence Dl— > D2) the fraction of colder gas 
in the mid-plane decreases and the fraction of hotter gas 
increases. From these two figures the following other 
features can also be highlighted: (a) the filling factors 
of the different gas components change strongly with the 
height above and below the mid- plane, \zd\- Typically, 
the fraction of colder gas decreases and that of the hotter 
gas increases with \zd\\ (b) run R2 is the only one where 
these filling factors remain nearly constant with \zd\, and 
it is also the only case that presents a non- negligible frac- 
tion of hot gas at low disk heights; (c) the ISM in the disk 
mid-plane of run Dl is composed only of the cold phase 
; (d) a multiphase ISM with a non-negligible fraction of 
cold gas at low heights exists only for run Rl. 

As can be seen in Figs. [12] and [T3l the only two runs 
that presents cold gas in the mid-plane are run Dl and 
Rl. The mid-plane of run Dl is dominated by cold gas 
because it has a high nsF value, a low stellar mass limit 
m-*,iim, and an inefficient stellar feedback. Run Rl, on 
the other hand, has a significant amount of cold gas in 
the mid-plane because in this run the UV background 
ionizing radiation was artificially reduced^ from z = 8 to 
present-day to its value at z — 8 and also because it has 
the largest (3, implying that some of the cold dense gas 
does not form stars, and therefore can remain as such. 

In Figure[14]we plot the vertical filling factor profiles of 
runs D3, D4, R2, and R3 at z = 0.43 in order to see how 
much these profiles could change with time. These are 
the cases that change most from z = 1 to z = 0.43. The 
general trend is that the fraction of colder gas increases in 
the mid-plane with time (specially in run R3), while the 
hotter gas 'moves' to higher vertical distances. In most 
of the runs, one observes that with time the warm gas 
fraction accumulates in an narrow extraplanar layer at 
~ 0.4 — 0.8 kpc that coexists with the hot gas. At larger 
heights, the fraction of hot gas already dominates. Run 

3 The formula for the ionizing flux, which is based on 
IHaardt & Madau I il99B) , is taken from lKravtsovl 1T999T) . For ref- 
erence, we note that the fluxes at redshifts z = 8, 2maxi 

and 1.0 

are 5.9 X 1CT 26 , 7.2 X 10" 22 , and 8.5 X 1CT 23 erg s" 1 cm" 2 Hz" 1 , 
respectively, where 2 max denotes the redshift at which the ionizing 
flux is maximum. 



D3 at z = 0.43 has at the mid-plane ~ 20% and ~ 80% of 
the volume occupied by cold and warm gas, respectively, 
while in run D4, ~ 100% is occupied by warm gas. The 
filling factor of the warm-hot gas in both runs peaks at a 
height ~ 700 pc, at ~ 75% and ~ 90%, respectively. The 
hot gas becomes the dominant component at a height 
~ 1 kpc. For the runs R3 and R2, the warm-hot gas is 
present in the mid-plane, with volume fractions of ~ 20% 
and ~ 40%, respectively. 




Fig. 12. — Fraction of volume filled with gas at different tem- 
peratures as a function of height (below and above the disk plane) 
for runs of series D at z = 1. All runs of the series D present no 
hot phase in the mid-plane of the disk (\z\ < 300 pc/h) while runs 
from D2 to D4 practically have no cold gas (T < 3000 K). On the 
other hand, run D4 is the only D model that has a fraction of hot 
gas that increases as we go away from the mid-plane of the disk. 
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Fig. 13. — As in Figure [12] but for runs of series R. As it was the 
case for runs D2 to D4, runs R2 to R4 have practically no cold gas 
whereas, as it happens with run D4, here runs Rl and R3 have a 
fraction of hot gas that increases as we go away from the mid-plane 
of the disk. On the other hand, run R2 is the only one in which we 
found a non-negligible fraction of hot gas in the mid-plane of the 
disk. 



12 



TABLE 3 
Global properties of models 



Model M vil M, f g a R e h f gal (R vir ) c f ga i(5R e ) d J/M(halo) J/M(baryons) 

(10 10 Mq) (10 10 M ) (kpc) (km s" 1 kpc) (km s" 1 kpc) 



2 = 1 



Dl 


4.69 


0.312 


0.412 


0.24 


0.16 





10 


96.0 


88.3 


D2 


4.69 


0.466 


0.070 


0.49 


0.12 





10 


203.9 


59.9 


D3 


4.87 


0.308 


0.446 


1.16 


0.13 





11 


208.8 


287.1 


D4 


4.87 


0.304 


0.210 


1.17 


0.10 





OS 


195.9 


143.3 


Rl 


4.79 


0.332 


0.133 


0.90 


0.13 





07 


188.1 


118.6 


R2 


4.99 


0.579 


0.076 


1.46 


0.14 





11 


205.2 


107.6 


R3 


4.81 


0.625 


0.043 


0.64 


0.15 





12 


218.9 


90.5 


Rl 


4.80 


0.524 


0.065 


0.37 


0.15 





10 


232.0 


51.3 










z 


= 0.43 










Dl 


5.40 


0.358 


0.396 


0.41 


0.17 





10 


143.9 


119.5 


D2 


5.54 


0.546 


0.060 


0.51 


0.13 





10 


218.1 


57.9 


D3 


5.84 


0.429 


0.369 


1.00 


0.13 





11 


199.6 


261.3 


Dl 


6.60 


0.309 


0.229 


1.11 


0.09 





06 


49.4 


149.7 


Rl 


5.64 


0.379 


0.161 


1.08 


0.13 





07 


240.8 


159.1 


R2 


6.21 


0.676 


0.075 


1.80 


0.13 





11 


141.2 


167.8 


R3 


5.49 


0.640 


0.043 


0.73 


0.13 





11 


203.7 


93.8 


Rl 


5.63 


0.571 


0.051 


0.50 


0.13 





10 


229.4 


49.9 



a Gas fraction as defined by f g = M gas /(M gaa + M„), where only gas cells with T g < 10 4 K are used. 
b Effective radius defined as the radius where half of the central galaxy stellar mass is contained. 
c Baryon fraction within the virial radius of the halo using all baryons, gas at all temperatures plus stars. 
d Baryon fraction within five effective radius using only stars plus gas with T g < 10 4 K. 
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Fig. 14. — Fraction of volume filled with gas at different temper- 
atures as a function of height (below and above the disk plane) for 
runs D3, D4, R2, and R3 at z = 0.43. These are the runs that 
most change their vertical filling factors since z = 1. 

3.5. The mass and angular momentum fractions 

In Table El M* and the gas fraction, f g , of the galaxies 
in the different runs are reported for z = 1 and 0.43. The 
gas fraction is defined as f g = M gas /(M gas + M*), where 
only gas cells with T < 10 4 K are used for the gas. The 
stellar and gas masses were measured within a sphere 
of radius 5R e . We see that for runs without efficient 
feedback, runs Dl and D2, f g strongly decreases (gas is 
consumed efficiently) as usf diminishes and, to a lesser 
degree, as m*jj m increases. Runs with a stochastic SF 
scheme do not show a clear trend of f g with tt-sf- From 
run R4 to R2, f g first decreases (runs R4— > R3) and then 
it increases (runs R3— > R2). On the other hand, run Rl 
is the run, within the R series of runs, with highest f g 



fraction at both z — 1 and 0.43. 

The shape of stellar particle mass distribution seems to 
influence the value of f g : for cases where this distribution 
extends to high-mass stellar particles, f g is smaller (from 
run D3 to D4 and partially from run Dl to D2). Note 
that in run D4 the feedback is so strong that a galactic 
wind is produced and as a result a signicant fraction of 
gas is ejected from the disk. Finally, from Table [3] it is 
inferred that runs with the random SF scheme are more 
efficient in transforming gas into stars than those with 
the deterministic scheme (the f g values of R runs are 
systematically smaller than those of runs D). However, 
it should be noted that as a consequence of the low value 
of C* and the SF scheme, gas is transformed locally into 
stars less efficiently in runs Dl to D3 than in the R runs. 

Columns (6) and (7) of Table [3] show the mass baryon 
fraction using all baryons (gas at all temperatures + 
stars) within i? V ir, and using only those baryons (gas 
with T g < 10 4 K + stars) within a sphere of 5i? e around 
the galaxy center, respectively. The latter case is related 
to what is called usually the galaxy baryon fraction. Sev- 
eral kind of estimates and observational inferences show 
that this fraction in present-day galaxies is much smaller 
than the universal value, equal to 0.15 for the cosmology 
used here. A radius of up to 5R e was used to compute 
f ga i to encompass the gaseous disk, which is typically 
larger than the stellar disk. 

The galaxy mass baryon fraction, f ga i{5R e ), does not 
depend on the value of ngp, as is seen from the runs 
where this parameter is changed systematically, (se- 
quences R4-» R3-» R2 and Dl-s- D2). From run D2 to 
D3, f ga i(J>Re) is almost the same, showing that switching 
off local cooling by 40 Myr, at least in this simulation, 
does not produce a net loss of baryons from the galaxy. 
However, from run D3 to D4, both the halo and galaxy 
baryon fractions decrease significantly. At z = 0.43 run 
D4 has the lowest values of these fractions among all the 
runs. Therefore, a bias toward massive particles m* dis- 
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tribution along with efficient feedback induces strong gas 
ejection from the galaxy and halo. 

According to Table |SJ the fraction of baryons within 
the halo but outside the galaxy is non-negligible for all 
the runs. The extreme cases are runs Dl, Rl and R4 
where 50% or more of the baryonic content of the cor- 
responding halo is outside the central galaxy. These 
baryons are in the form of shock-heated gas or gas heated 
and expelled from the central galaxy; a small fraction is 
in stars in satellite galaxies. The halo with the smallest 
baryon fraction corresponds to run D4 at z = 0.43, with 
~ 60% of the universal fraction. This run is the one with 
the most efficient galaxy outflows. 

Columns (8) and (9) of Table [3] present the total spe- 
cific angular momentum (AM) for the dark matter halo 
and for the sum of the stellar and cold gas components 
of the central galaxy (within a sphere of 5R e radius), 
respectively. The specific AM is given in units of km 
s _1 kpc. The total specific AM measured at a given 
time is sensitive to the current halo and galaxy assembly 
events, nevertheless, with a couple of exceptions, values 
are around 200 km s _1 kpc, a value similar to those found 
in the adiabatic run (see Table [2J . In this latter case the 
halo is almost not 'perturbed' by the baryonic processes. 
We found that the specific AM of the DM halo is typcally, 
as expected, greater than the corresponding specific AM 
of the baryonic component. The effects of complex bary- 
onic processes in the different runs produce variations in 
the halo mass assembly history as well as in the prop- 
erties of the surviving subhalos. Such variations are on 
the basis of the different values of the halo's specific AM 
measured in our different runs, at epochs where the halo 
is yet in a relatively active phase of mass aggregation. 

In spite of the transient nature of the specific AM mea- 
sures mentioned above, some general results can be high- 
lighted. We see, as expected, that in most of the runs 
the halo's specific AM is larger than the galaxy's spe- 
cific AM. Nevertheless, the galaxy's specific AM for those 
runs with strong feedback is relatively large and in some 
cases it is even larger than the corresponding halo's spe- 
cific AM (see Table H run D3). The specific AM in most 
of the runs remains roughly constant or slightly increases 
from z = 1 to z = 0.43. The effect of feedback is clearly 
in the direction of producing galaxies with larger specific 
AM (see the change from run D2 to D3). We also see 
that the galaxy's specific AM increases as nsF decreases 
(sequence R4— > R3— > R2), but this seems to be related to 
strength of the feedback, it is stronger in low gas density. 

4. DISCUSSION 

4.1. General trends and SF/ISM properties 

In the sequences of runs D2 — )• Dl (deterministic SF 
model with weak feedback) and R2— > R3— > R4, (ran- 
dom SF model with efficient feedback) the SF density 
threshold hsf is increased, keeping the other parameters 
similar. Our simulations show that as uqf is increased, 
the overall stellar galaxy has a smaller effective radius R e 
and the circular velocity profile is more peaked. For lower 
values of nsF, the peak in the SFH tends to be slightly 
shifted to later epochs and the later (nearly exponential) 
decline in SFR is more gradual. Related to this trend, at 
smaller nsF, the stellar populations become younger on 
average, the stellar structure becomes thicker (see also 



iSaitoh et al.l [20 081 . and a larger fraction of stars exists 
in a "puffier", more spheroidal component. 

Regarding the ISM properties of the simulated galax- 
ies, their PDFs of the density and temperature evolve to 
broader distributions, with more gas in the low-density, 
and high-temperature regime in the disk, and a smaller 
maximum value for the gas density, as tisf is smaller. 
The gas filling factors above and below the mid-plane 
also show that for smaller rtgF, the cold gas close to the 
mid-plane become less abundant. We note that in gen- 
eral is not easy to produce cold, high-density gas in the 
simulations. At z = 1.0, only the run with the high- 
est tt-sf (= 50 cm -3 ) and weak feedback (run Dl) shows 
a significant amount of this gas in the mid-plane. At 
z = 0.43, cold and high-density gas is also present in the 
mid-plane of run D3. On the other hand, a hot gas phase 
in the disk mid-plane can only apparently be produced 
by lowering nsF- 

It is not clear what a realistic value for n$F should be, 
but we argued in §2.2 that for a cell of 150 pc, close to 
the 218 pc of our nominal resolution at z = 0, this value 
should be around 5 cm~ 3 . Moreover, the results of our 
experiments seem to support this number: a value much 
higher than this produces a galaxy that is too compact 
(see, for example, R2 and R3 versus R4 or D2 versus 
Dl, albeit with weak feedbak), while a much lower value 
produces a system that is closer to being a spheroidal 
galaxy, with a very small fraction of gas, rather than a 
disk (see Table [3]). However, a small value of usf (0.1 
cm~ 3 ) ma y produce reasonable res ults as shown, for ex- 
ample, by IGovernato et aLl (|2007D . We speculate that 
this may be an effect of the different technique (AMR 
versus SPH) and subgrid physics (feedback) used to sim- 
u late galaxy fo r matio n. 

ISaitoh et al.1 (|2008ft concluded that values of nsF as 
high as ~ 100 cm~ 3 should be used in high-resolution 
simulations (~ 20 pc) in order to reproduce the complex 
structure of the gas disk. They also concluded that us- 
ing high values for usf makes the model results fairly 
insensitive to the SF prescription. In our experiments, 
high values of nsp produce systematically too concen- 
trated galaxies with the hot gas phase being absent from 
the mid-plane. Yet, results in ISaitoh et al.l (|2008f) are 
actually in agreement with our suggested strategy for 
choosing nsF on the basis of the gas column density. For 
a 20-10 pc resolution, our simple argument would indeed 
suggest values of n$F of 50-100 cm -3 

A second parameter that enters in the SF prescription 
is the stellar particle mass limit TO*ji m , defined in §2.2. 
As was mentioned in that section this is not an abso- 
lute lower limit. High values of to*, for a given cell size, 
mean that more stellar and SN energy is dumped into 
the cell. As a result, in models with cooling turned off in 
the SF regions, the feedback-driven ejection of gas from 
the disk increases and the SFH of the modeled galaxy 
becomes more bursty (compare runs D3 and D4). On 
the other hand, the structural/dynamical properties of 
the modeled galaxies does not seem to depend strongly 
on the stellar particle mass distribution. 

Regarding the efficiency of stellar feedback, we find 
that it has a strong effect on the structural/dynamical 
properties of the galaxy as well as on their SFH and ISM 
properties, in agreement with recent numerical works 
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Fig. 15. — Edge-on view of simulation R3 at z = 1. Left panel: Distribution of the gas temperature. Right panel: Distribution of the z 
velocity component. The slices are 36 comoving kpc on a side (length of solid line) and 219 comoving pc thick. It can be seen that gas 
velocities can clearly exceed 200 km s" 1 . 



c.f. iGovernato et al.l 12004 120071: iScannapieco et al 



255^ , I2008t iCeverino fe Klypin I l2009t iZavala et al 
IDalla Vecchia fc Schave 2008; Stin son et al.H2OO90 



2008; 



We have experimented by switching off cooling (c.f . 



Gerritsen fc Ickel [l997t iThacker fc Couchman I T2000t 



of runs Dl, D3, D4, and R3 at z = 1. Here we have used 
all cells inside i? V ir- Unlike runs Dl and D3, runs D4 and 
R3 exhibit high \v z \ wings, corresponding to a galactic 
wind. 



Stinson et al. 2006) by 40 Myr in order to allow for 



proper expansion of the heated regions. This time is 
roughly the time a pressure-driven super-shell will at- 
tain > 150 pc size and the time a star of eight solar 
masses live. When this feedback is allowed in our low- 
mass galaxy, the stellar effective radius increases by a fac- 
tor of ~ 2, the galaxy becomes less concentrated, and its 
V C (R) profile becomes flatter than in the weak- feedback 
simulation (c.f. runs D2 and D3). A stronger feedback 
helps also to keep an episodic SFH. A complex multi- 
phase structure in the disk is favored by the action of 
feedback but if it is too strong, the low-density, warm 
and hot gas can be completely blown away from the disk, 
as it seems to be the case in run D4. 

The warm and hot gas seen above the disk plane can 
be produced by two mechanisms: (1) shocks of the in- 
falling gas, and/or (2) gas heated and ejected from the 
mid-plane by supernovae and stellar winds. Which is the 
dominant process depends on the infalling history of the 
gas and on the efficacy of the stellar feedback. For ex- 
ample, in run Dl (weak feedback) it is very unlikely that 
the hot gas seen at intermediate altitudes (~ 1 — 3 kpc) is 
produced by the second mechanism. On the other hand, 
the increasing amount of hot gas seen at very high al- 
titudes (as high as > 5 kpc) in runs D4, R3, and R4 
at z = 1 appears to come from the second mechanism. 
This is confirmed visually and by the analysis of the gas 
v z distribution, where v z is the gas velocity component 
perpendicular to the plane of the disk. Figure [15] shows 
run R3 at z = 1. The projected spatial distribution of 
the gas temperature and of the z velocity component arc 
shown in left and right panels, respectively. Figure 1161 
on the other hand, shows the v, normalized distribution 




-200 200 
V, (km/s) 

Fig. 16. — Distribution function, weighted by volume, of the 
vertical velocity component of the gas for the analyzed galaxy at 
z = 1 for runs Dl, D3, D4, and R3 . We have used all cells within 
R v ir- As was already seen in Fig. ll5l runs D4 and R3 exhibit outflow 
velocities that go up to 400 km s _1 . 



4.1.1. Self -regulated SF and the multi-phase ISM 

Several of the results described above can be under- 
stood as a consequence of the efficiency by which the 
thermal energy is injected into the ISM and the ability of 
the medium to self-regulate its SF. The R2-R4 random- 
SF sequence (cooling turned off after SF events) shows 
that the reduction in usf causes: (a) a flatter circular 
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velocity profile; (b) a stronger spheroidal stellar compo- 
nent; (c) a thicker disk; and when run R2 is compared 
with R3 or R4, (d) a delay of the maximum of SFR. 
The deterministic-SF sequence shows that, at fixed usf 
(runs D2, D3, and D4), a greater ability of SF to heat 
the gas causes a more efficient self-regulation of SF, pre- 
venting the gas from going into stars too rapidly. This 
gives flatter V C (R) profile and allows the gas to settle into 
an extended disk configuration. However, when feed- 
back is strong enough to produce galactic superwinds, 
the self-regulation can be interrumped and SF proceeds 
in a bursty regime. Indeed, run D4 is the one where 
ejecting feedback is most efficient and it is also the one 
with the most bursty SFH. 

A multi-phase structure of the galactic ISM is not easy 
to attain in cos mological numerical s imulations (see for a 
discu ssion, e.g.. iStinson et al]l2006HCeverino fc Klypinl 
2009). For example, we see in run Dl that only the cold 
phase is present in the mid-plane of the disk while in 
run D2 the hot phase is absent (see Figure IT2j) . These 
two runs have an inefficient feedback, but even in those 
for which the cooling is switched off for 40 Myr after 
each SF event (strong feeback, runs D3 and D4), the 
hot phase is also absent. The only run in which a a 
multi-phase medium is observed at the analyzed epochs 
is R2. Its density and temperature PDFs show wide and 
nearly bimodal distributions, with different gas phases, 
from cold-warm dense gas to hot low-density gas, coex- 
isting within the galaxy. The increase of the SFR with 
time of run R2 is more gradual and not too bursty, which 
suggests a more uniformly self-regulated SF process with 
only small episodic phases. Finally, this run takes the 
longest time to assemble the disk. 

In summary, the results of our experiments show that 
a multi-phase ISM in low-mass galaxies is not a straight- 
forward product of feedback and resolution. While the 
former certainly helps, in low-mass systems it also pro- 
duces outflows that leave only cold/ warm high-denstiy 
gas in the galaxy and induce a bursty SFH. 

Finally, in our simulations the Esfr - ^ relation tends 
to be significantly steeper than the empirical determina- 
tions for normal disk galaxies, which are typically ex- 
plained by the quiescent, self-regulated nature of the 
SF process (e.g..lKennicutt 111 9981: Hernandez et al.ll2001t 
Marti n fe Kennicuttll200lI TFo r reference, we plot in Fig- 
ure [7| the fit to observations bv lKennicutt I (|1998f ) for av- 
erages values of normal disk galaxies, their centers, and 
starbust galaxies (solid line) . This correlation is roughly 
maintained for azimuthally averaged elements along a 
given disk (|Martin fe Kennicuttil200ll ). As can be seen, 
most of our runs at both z = 1 and z = 0.43 show a 
steeper local correlation than the observational fit of nor- 
mal galaxies. Overall, most of the simulation points lie 
around the observational fit. The points from run Dl, 
which is the run that suffers most from the overcooling 
problem, are the most distant from the fit. It is inter- 
esting to note that among LSB galaxies and along their 
azimuthally averaged elements, the £sfr-£ s correlation 
has been found to be significantly st eeper than the given 
by Kenniccutt (jWvder et al. 1 120091 ). LSB galaxies are 
typically of low mass. 

4.1.2. The simulation with run-away stars (run Rl) 



The subgrid physics used in run Rl differs from the 
rest of simulations explored herein. In run Rl the radia- 
tive cooling is not switched off but a signicant fraction 
of stars are modeled as runaway stars able to migrate 
to lower gas density regions. The SN explosions in the 
lower density gas, where the cooling time is long, are 
more efficient. This model has been succesful in repro- 
duc ing many of the properties of Milky Way-sized g alax- 
ies HCeverino fc Klypin Il200l iCeverino et al. 1120091) . As 
shown here, it works reasonably well also for a low-mass 
galaxy. It successfully reproduces a roughly flat V C (R) 
profile and a multiphase ISM at the galaxy mid-plane. 
The model also has the lowest galaxy mass baryon frac- 
tion and relatively high angular momentum content. On 
the other hand, its SFR was relatively too efficient at 
early epochs (the first 1-2 Gyr). As a result, at later 
times, z < 1, the SFR was lower than other 'feedback- 
efficient' models, and the model ends up with a lower 
fraction of mass in gas (see, for example, Rl versus D3). 
These effects are probably caused by the high value used 
for the star formation efficiency, C*. Overall, run Rl ex- 
hibits similar properties to those of run D3 and D4 as a 
result of its ability to heat the gas efficiently through the 
runaway stars. 

4.2. Galaxy assembly in low-mass ACDM halos 

Our simulations suggest that the effects of feedback are 
more relevant for the galaxy morphology than in massive 
systems. The main effects are galaxy outflows and a trend 
to form a significant extended stellar spheroid with an ap- 
proximately exponential surface density profile. What is 
the origin of this stellar halo? Since the dark matter 
halo under study evolves in a quiet fashion, at least from 
z ~ 2 to z ~ 0.4, without experimenting major merg- 
ers, a significant fraction of the stellar particles of this 
extended spheroid were born in situ and not by mergers 
as it is expected for large galaxies. This is confirmed by 
our analysis of tracing back in time the position of the 
stellar particl e s in o rder to infer their origin. Recently, 
IStinson et al.1 (|2009( ) have also reported the formation of 
old extended exponential stellar halos in their simula- 
tions of dwarf-like galaxies in isolated (without mergers) 
halos. As these authors review, observations show that 
dwarf galaxies are characterized by extended exponen- 
tial old/intermedium age stellar halos. A recent study 
of LSB galaxies with the stacked image technique has 
also shown the existence of red faint stellar halos around 
LSB galaxies (jBergvall et al. 1 12009). galaxies that are 
supposed to form typically in isolated environments. It 
is worth noticing that the R runs (those with a stochas- 
tic SF model), with the exception of R2, produce more 
massive and extended stellar halos than the D runs (those 
with deterministic SF model), as seen in Fig. [3] 

The parameters of the SF/feedback prescriptions used 
here were shown to affect significantly the evolution and 
properties of the modeled low-mass galaxy. The change 
of these parameters affects in a complex and non-linear 
way the evolution and properties of the system. We note 
that even the evolution of the dark matter mass and 
angular momentum distributions is different for differ- 
ent combinations of the SF/feedback parameters, con- 
firming the existence of a tight evolutionary intercon- 
nection betw een baryons and dark matter (see for recent 
results, e.g.. I Weinberg et al. Il2008t lAbadi et al. 1 120091: 
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Fig. 17. — Edge-on view of model D4 at z = 0.43. Left panel: Distribution of the gas temperature. Right panel: Distribution of the 
velocity component perpendicular to the disk. Slices are 25 comoving kpc on a side (length of solid line) a nd 2 19 comoving pc thick. It can 
be seen that the gaseous disk is warped and more homogeneously distributed than that shown in Figure [T5l Gas velocities are no longer 
structured in shells and velocities exceed 300 km s _1 far from the disk. 

the rotational velocity, of the gas at T < 10 4 K, increases 
as observed in low -mass galaxies. Similar situations have 
been discussed bv lValenzuela et al"~l (|2007l ) as a result of 
pressure gradients and non-circular motions. The mis- 
match between circular and rotation velocity deserves 
further exploration in future papers exploring also effects 
of resolution. 

Also, our simulated galaxies show baryonic mass frac- 
tions that are too high compared with observational in- 
ferences. For example, in the simulation where galaxy 
outflows are most effective (run D4), f ga i(5R e ) — 0.06 
(at z — 0.43) is only a factor 2.5 lower than the uni- 
versal baryon fraction while the observational infer- 
ences suggest that present-day galaxies of low masses 
have baryon fractions more than a factor of 5 smaller 
than /{,. Our results confirm that 'producing a stellar 
feedback that is sufficiently efficient as to eject consider- 
able amounts of gas from even low-mass galax ies is not 
an easy task (see also lDubois &: Tevssier ll2008l and ref- 
erences therein). 

Even more, the problem of producing strong enough 
outflows is probably in opposition with other potential 
vexing conflict. An increasing amount of local and high- 
rcdshift observational evidence shows that the less mas- 
sive the galaxies, the higher their current specific SFRs 
(SSFR = SFR/M,) on average at all rcdshifts back to 
z ~ 1 — 2, a phenomen on dubbed as 'downsizing in SSFR' 
(|Fontanot et al.l [2009(1 . Moreover, the observational in- 
ferences show that the problem is not that massive disk 
galaxies have low SSFRs but that lo w-mass galaxies have 
too high SSFRs (see for a discussion lFirmani et al. II2009L 
and references therein). Such SSFRs imply that low- 
mass galaxies should delay the formation of stars and 
do it with a nearly constant or even increasing SFRs 
(|Noeske et al. 112007ft . On the other hand, if f gal de- 
creases as the galaxy mass decreases, then obtaining late, 
high SSFRs in these galaxies becomes even more diffi- 



iPedrosa et alJl2009| ) . In particular, the internal distribu- 
tion of the angular momentum evolves in our simulations, 
as can be seen, in particular, by changes in the dark halo 
and baryonic galaxy specific angular momentum ratio. 
We suggest that our finding that the baryon specific an- 
gular momentum can be larger than the one in the dark 
halo, deserves further future exploration. 

In order to obtain a not too concentrated low-mass 
galaxy, with a nearly flat rotation curve, multi-phase 
ISM, and an extended and episodic SFH, an efficient 
stellar feedback is necessary (in our case, switching off 
radiative cooling in the star forming region by 40 Myrs) . 
However, the SF parameters are also important; for the 
resolution we have attained in our simulations (nomi- 
nally 218 comoving pc), effective values for the density 
threshold usf ~ 1 — 10cm~ 3 are recommended. We say 
'effective' because in the stochastic (random) SF scheme 
the reported tisf is just the minimum value, larger val- 
ues of nsF are actually chosen, and SF proceeds with a 
probability distribution given by eq. ([T|). As in previ- 
ous related numerical treatments of the galaxy problem, 
we conclude that sub-grid physics needs to be carefully 
addressed (tuned) in order to obtain reasonably realistic 
galaxy properties. 

4.3. Shortcomings 

Despite that several of the runs presented here predict 
realistic properties^, we note that the circular velocity 
curves arc decreasing for all the cases studied here, while 
in observed low-mass galaxies the trend is to have an in- 
creasing rotation curve up to the last measured radii (c.f. 
IPersic et all 119961 but see Swaters et al. 2009). How- 
ever, the average rotational velocity profile does not trace 
necessarily the circular velocity and at least in model R2 

4 The aim of the present work was not to produce a realistic 
galaxy. 



17 



cult . As concluded in [Firmani et al. I (|2009l see also 
Somerville et al. 2008 and Fontanot et al. 2009) both 
problems are difficult to solve in the context of ACDM- 
based models. 

The cosmological numerical simulations presented here 
confirm the sharpness of both potential problems. The 
baryon fractions of the simulated low-mass galaxies are 
too high and their SSFRs are in all cases much smaller 
than the average values of the observational inferences at 
redshifts ~ 0.4 - 1. The run with the highest SSFRs is 
D3: SSFR = 0.17, 0.06 Gyr" 1 at z = 1, 0.43, respec- 
tively. Observations of large samples of late-type galax- 
ies at high redshifts report values ~ 5 ti mes higher for the 
corresponding z's and stellar mass es dBell et al. I [20071: 



2003J iMac Low fc Klessenl 120041: iRosas-Guevara et all 
2010 ; IVazquez-Semadeni et al.l I2010D . The magnetic 
field is expected to h a ve a si m ilar effect (Ostrikcr et al. 
1999[ iHeitsch et all 120011 : IVazquez-Semadeni et~ 
2005). Note that, in a sense, this is what we are 



Zheng et al.Tl2007t iNoeske et al~l l2007t iDamen et al.1 



2009L see also Fig. 1 in Firmani et al. 2009). Most of the 



other runs have SSFRs much smaller than those of run 
D3. Our simulations also show that if f ga i is reduced 
by strong feedback-driven outflows, then the late-time 
SSFR decreases. Run D4 indeed has fewer baryons ow- 
ing to the strong outflows, but also a much lower SSFR 
than run D3. A related problem to the low SSFRs is that 
the gas fractions of the simulated galaxies (for example 
at z = 0.43,) are in most cases too low (see Table |3|). 
Among runs R2-R4, the maximum value of f g is 0.075 
for run R2; among runs D2- D4, as expected by its SSFR 
history, run D3 has the largest value, f g = 0.37, which 
is still marginal compared with the large gas fractions 
of observed present-day low-mass galaxies. They should 
have large gas reservoirs to produce the high SSFRs in- 
fered from observations. 

If the baryon fraction and SSFR problems for low-mass 
galaxies are confirmed, models for SF/feedback very dif- 
ferent to the one applied here — and in other works — 
should be proposed to agree with observations, preserv- 
ing the underlying ACDM cosmogony. It is possible that 
significant changes to the galaxy assembly process can be 
introduced alternatively by improving the physics of the 
poorly understood circumgalactic and intergalactic me- 
dia. Finally, we should say that the current observational 
inferences, specially at high redshifts, are still plagued by 
large uncertainties (iConrov et al.l 120 09a: Co nrov et al. I 
I2009bt iFhmani et al. 1120090 . 

A comparison of simulation outcomes at z = with 
the more accurate observations of present-day galaxies 
is desirable in order to asses the severity of the problems 
mentioned above. We plan to do this in a forthcoming 
paper. 

Finally, it should be mentioned that the sub-grid 
prescriptions in our simulations are certainly not 
complete, as several physical ingredients that could 
be important for the evolution of real galaxies are 
missing. In particular, our simulations do not include 
the magnetic field nor cosmic ray pressure, and they 
do not resolve the whole plethora of phenomena that 
occur at the scales of molecular clouds, including 
turbulence and other forms of stellar feedback additional 
to supernova explosions, such as bipolar outflows and 
the expansion of HII around massive stars. All of 
these agents contribute to drive the turbulence at 
scales that are unresolved in our simulations, whose 
effect may be to r educe the star formation efficiency 
in the simulations dKlessen. Heitsch fc Mac Low! [2 000: 
Vazquez- Semadcni, Ballcstcros-Parcdes, & Klessen 



m a sense, 

accomplishing in our simulations by turning off the local 
cooling after a star formation event, since the hotter gas 
is less prone to forming stars. However, future studies 
will need to quantify these effects so that they can be 
incorporated in a suitable sub-grid model applicable for 
galaxy-formation simulations. 

5. CONCLUSIONS 

By means of high-resolution AMR hydrodynamic sim- 
ulations (the hydrodynamics + N-body ART code), we 
have explored the formation and evolution of a galaxy 
within a low mass halo (~ 7 x 10 10 h~ 1 MQ at z = 
0) assembled in the context of the ACDM cosmogony. 
The DM particle mass of the least massive species is 
5.3 x 10 5 /i _1 M©. Our main aim was to explore the 
formation of a low-mass galaxy in cosmological simu- 
lations with the same initial conditions but varying the 
most relevant parameters of the SF and feedback recipes. 
The analysis presented herein was mainly made at two 
epochs, z = 1.0 and 0.43, a time period of 3.2 Gyr for 
the used cosmology. From this exploratory study -not 
exaustive in any sense- we remark the following system- 
atic results: 

- The SF density threshold tt-sf is a relevant parame- 
ter in the simulations, which affects the main structural 
and dynamical properties of the stellar and gaseous com- 
ponents of the low-mass galaxies, as well as their SFHs 
and general ISM properties. As n$F is decreased, (a) the 
stellar galaxy size (R e ) increases; (b) the V C (R) profile 
becomes less peaked; (c) the peak in the SFH is delayed 
and the late (nearly exponential) decline in SFR becomes 
more gradual; (d) the stellar populations are younger on 
average, with a thicker disk and a larger fraction of stars 
in a "puffier" , spheroidal component; (e) the low-density 
and high-temperature gas in the disk increases. For a cell 
of 150 pc, close to our nominal resolution at z = Q, we 
estimate from typical ISM properties that tisf should be 
~ 5 cm~ 3 . 

- The stellar particle mass limit, m» ji m , affects the effi- 
ciency of SN energy dumping into the cell. For high val- 
ues of m»ji m , the feedback-driven disk outflows increase 
and the SFH of the simulated galaxy becomes burstier. 

- An effective stellar feedback (generated by turning 
off locally the cooling by 40 Myr after SF) has a strong 
effect on our studied low-mass galaxy: the stellar radius 
R e increases by a factor of ~ 2, the galaxy becomes less 
concentrated, and V C (R) becomes flatter. A strong feed- 
back produces also an episodic SFH. Its action also favors 
the development of a complex multi-phase disk ISM and 
a self-regulated SF process, but if it is too strong, the 
low-density, warm/hot gas can be completely blown out 
from the disk. 

- For all of our runs, both at z = 1 and 0.43, we find 
that the relation between Ssfr and S s is steeper than 
empirical determinations for normal disk galaxies, but 
in agreement with recent determinations for LSB galax- 
ies, which are typically of low mass. Nevertheless the 
effective stellar surface densities, Em, of our models are 
instead in the high-density side of the observed distribu- 
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tions at z ~ 0.4 and 1. 

- In all runs, the stellar galaxy has a significant ex- 
tended spheroidal component (between ~ 30 and 60% of 
M* at z = 0.43) with a nearly exponential surface den- 
sity profile. We have seen that an important fraction of 
the stars in this halo were formed in situ, while a minor 
fraction may come from minor merger events. 

- The variation of the parameters of the SF/feedback 
prescriptions affects in a complex and non-linear way 
the evolution of the galaxy system. Even the evolution 
of the dark mass and angular momentum distributions 
change as these parameters are changed, confirming the 
tight interconnection between baryons and dark mat- 
ter. For example, the internal distribution of AM evolves 
in our simulations and this evolution is sensitive to the 
SF/feedback physics. 

- In all runs of the low-mass galaxy studied here, the 
baryon fractions are higher, the gas fractions are lower, 
and the SSFRs are much smaller than the corresponding 
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